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Abstract 

The upsilon distribution, the sum of independent chi random variates 
and a normal, is introduced. As a special case, the upsilon distribution in¬ 
cludes Lecoutre’s lambda prime distribution. [22] The upsilon distribution 
finds application in Frequentist inference on the Sharpe ratio, including 
hypothesis tests on independent samples, confidence intervals, and pre¬ 
diction intervals, as well as their Bayesian counterparts. These tests are 
extended to the case of factor models of returns. 


1 Introduction 

For k > 1, given a fc-vector 
pi, 1/2, ■ • -^VkV, if 

y 

where Z ^ Af {0, 1) independently of Xi ^ which are independent, then 

we say y follows an upsilon distribution^ with coefficient t = ■ ■ ■ ,tk] , and 

degrees of freedom v> = pi, 1 / 2 ,..., z/fe] . Let this be written as y ^ T (t, i/). 
The summands are reminiscent of Nakagami variates, but the latter may only 
take positive weights by definition, whereas the elements of t are unrestricted. 

[14] 

Lecoutre’s lambda prime distribution appears as a specific case of the upsilon, 
with k = \. [15, 22] It should be stressed that the use of the lambda prime for 
inference with the t-statistic was first described independently by Akahira and 
Lecoutre and the results here should be viewed as extensions of that work. 
[1, 2, 16] Lecoutre does not consider the Sharpe ratio, which is simply the t 
statistic rescaled, nor does he consider the multiple independent samples case, 
which is first described here. 

* spavQalumni.emu.edu 

have chosen T since it vaguely looks like a ‘y’, half way between the ‘x’ of the chi 
summands and the ‘z’ of the normal. The ‘t’ is used for the coefficient since in practice the 
elements of this will typically be t-statistics. 


[ti, ^ 2 , ■ • ■, f/c]^ and fc-vector of positive reals 
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1.1 Computing the upsilon distribution 

It is customary, when introducing a new probability distribution, to quote its 
density, and describe cases in which it arises in real settings, how to interpret 
its parameters, etc. The upsilon is somewhat idiosyncratic in that it is not a 
distribution of primary interest, however, and only useful as a kind of ‘dual’ 
to the distribution of the Sharpe ratio. [27] As such, even though one could 
write down the density, as the convolution of the densities of the normal and 
the chi distributions, doing so would not much help in practical computation of 
the density. Moreover, inference on the parameters t is not typically required. 

However, computation of the cumulative distribution and the quantile of the 
upsilon are of interest in performing inference on the Sharpe ratio. These are 
complicated to write down: the CDF is the integral of the ugly convolution 
alluded to above, and the quantile requires inverting that integral. It would 
not seem a fruitful pursuit to find highly accurate methods of computing these 
functions, since the tests described herein for inference on the Sharpe ratio are 
only exact in the case of Gaussian asset returns, which is a poor approximation 
for most real assets. Instead, methods for approximating the distribution and 
quantile function are likely to be sufficient for most users. 

Approximating these functions via the Edgeworth and Cornish-Fisher ex¬ 
pansions is fairly simple, and relies only on knowing the cumulants of the up¬ 
silon distribution. [3, 17, 11] These are readily available since the cumulants 
commute with addition for the case of independent random variables. This 
means one only needs to know the cumulants of the normal distribution^ and 
the chi distribution. Users of the R language have access to the approximate 
density, distribution, and quantile functions via the sadists package available 
from GRAN. [23, 21] 

For concreteness, let ki be the raw cumulant of the upsilon distribution 
with coefficient [ti,t 2 , ■ ■ ■ ,tk]^, and degrees of freedom [-^ 1 ,^ 2 , ■ ■ ■ ■ Note 

that Ki is the expected value, and K 2 is the variance of this distribution. Then 


Ki 


K2 


= V2 E 


i<j<^ 


k3= E 



ti 

= 1+ E w 


— 

i<i<fc ^ 


( 2 ) 


In practice, it is easier to first compute the raw moments of the chi distri¬ 
bution with Vj degrees of freedom, since the raw moment is simply 



^The only distribution with a finite number of non-zero cumulants. 
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One then can translate the raw moments of the chi variables to raw cumulants. 
The cumulants of the upsilon can then be computed using properties of the 
cumulants (invariance, additivity and homogeneity). 


2 The Sharpe ratio 

The Sharpe ratio is the most commonly used metric of the historical perfor¬ 
mance of financial assets, t is defined as 


p A - 

C =df —:— 

a 


(3) 


where A is the historical, or sample, mean return of the mutual fund, a is the 
sample standard deviation of returns, and rg is some fixed risk-free or disastrous 
rate of return. Under Sharpe’s original definition, rg was equal to zero. [26] One 
typically uses the vanilla sample mean and the Bessel-corrected sample standard 
deviation in computing the ratio, 


M =df 


^ ^ (7 —df 

1 <t<n 


IE 


l<i<n 


(xt - A) 


n — 1 


(4) 


where Xt are the returns of the asset at time t. Herein, Sharpe ratio will refer to 
this quantity, computed from sample statistics, whereas signal-noise ratio will 
refer to the analagously defined population parameter. 


C =df 


M - ^0 
a 


(5) 


In general, hats will be placed over quantities to denote population estimates. 

Interpretation of the Sharpe ratio relies on a very simple model of returns: 
identically distributed, homoskedastic, and unconditional on any state variables. 
A more realistic model of real assets allows attribution of returns to contempo¬ 
raneous observable factors. In the general case one attributes the returns of the 
asset in question as the linear combination of I factors, one of which is typically 
the constant one: 

i-i 

Xt = Aol + Pifi,t + €t, (6) 

i 


where fi^t is the value of some A*' ‘factor‘ at time t, and the innovations, e, are 
assumed to be zero mean, and have standard deviation a. Here we have forced 
the zeroth factor to be the constant one, /g_t = 1. [25] 

Given n observations, let F be the n x I matrix whose rows are the observa¬ 
tions of the factors (including a column that is the constant 1), and let x be the 
n length column vector of returns; then the multiple linear regression estimates 
are 


/3=df(FTF) 


CT =df 


\ 


(a; - F,3) (a: - F^) 


.-I 


(7) 


We can then define a Sharpe ratio for factor models as follows: let v be some 
non-zero vector, and let rg be some risk-free, or disastrous, rate of return. Then 
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define 


Cg —df 


$ V -To 


( 8 ) 


Typically in the study of equity strategy returns, one chooses v = eg, the vector 
of all zeros except a one corresponding to the intercept term. 

There are numerous candidates for the factors, and their choice should de¬ 
pend on the return series being modeled. For example, one would choose dif¬ 
ferent factors when modeling the returns of a single company versus those of 
a broad-market mutual fund versus those of a market-neutral hedge fund, etc. 
Moreover, the choice of factors might depend on the type of analysis being per¬ 
formed. For example, one might be trying to ‘explain away’ the returns of one 
investment as the returns of another investment (presumably one with smaller 
fees) plus noise. Alternatively, one might be trying to establish that a given 
investment has idiosyncratic ‘alpha’ (ie., /3q) without significant exposure to 
other factors, either because those other factors are some kind of benchmark, 
or because one believes they have zero expectation in the future. [6, 13] 


2.1 The Sharpe ratio under Gaussian returns 


While normality of returns is a terrible model for most market instruments [7], 
it is a terribly convenient model. We will adhere to this terrible model and 
assume returns are unconditionally i.i.d. Gaussian, i.e., Xt ~ A/'(/i,cr^). The 
equivalent assumption for the factor model. Equation 6, is that the innovations 
are i.i.d. Gaussian, i.e., Ct ~ W (O, tr^). 

Under Gaussian returns, the Sharpe ratio follows a t distribution up to 
scaling. That is 


y/nl 


Z + V^C 


(9) 


where Z is a standard normal independent of the chi-square which has 

n — 1 degrees of freedom. The Sharpe ratio under the factor model, defined in 
Equation 8, has a considerably more complicated distribution; 




z 




- 1/2 


xl-iKn-l) 


( 10 ) 


where Cg = ^ is the population analogue of (g , and again Z is a standard 

normal independent of Xn-/’ ^ chi-square with n — I degrees of freedom. For 
the exact results described in this paper, one must make further unrealistic 
assumptions regarding the factor model, namely that the factor returns are 
essentially deterministic. 


3 Prequentist inference on the Sharpe ratio 

Assuming normal returns, the connection between the Sharpe ratio and the up- 
silon distribution becomes apparent when one rearranges the terms. Conditional 
on observing 

y/n( = Z + xl-i/in - 1). 
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(Since Z is an unobserved standard normal random variable, we will often flip 
its sign without renaming it.) Thus conditional on v^C takes an upsilon 


distribution with k = 1, t = y/n( and degrees of freedom v = [n — 1]. The 


same rearrangement can be performed for multiple independent samples, leading 
to tests of linear combinations of the population parameters. 

3.1 Hypothesis tests 

The classical one-sample test for the signal-noise ratio could be interpreted in 
terms of the upsilon, although it is far more typical to rely on the t-distribution 
instead. The typical way to test 


Ho:C = Co versus ili ; C > Co 


at the a level is to reject if the statistic t = y/n( is greater than ti-a — 1), 


the 1 — a quantile of the non-central t-distribution with n — 1 degrees of freedom 
and non-centrality parameter 6o = \/nCo- Alternatively one can perform this 



test by rejecting at the a level if v^Co < Ta ~ ij : the a quantile of 

the upsilon with t = i/nC u = [n — V\. 


For the one sample case, use of the upsilon is superfluous^, so consider the 
case of testing on multiple independent samples. A two-sample test for equality 
of signal-noise ratio, given independent observations, appears, at hrst glance, 
to be related to the Behrens-Fisher problem, which has no known solution. [5] 
However, this hypothesis can be tested exactly. For i = 1,2, given rii i.i.d. 
draws from Gaussian returns from two assets with signal-noise ratios Q y to test 


^^0 : Cl = C2 versus i?i ; Ci > C2, 


compute the sample Sharpe ratios, Cz- Again, rearrange to get 



where the Zi ~ A/'(0,1) independently and independent of the chi-square ran¬ 
dom variables Xi ^ (^i ~ l)i which are independent. 

Then under the null, 



where Z ^ Af {0, 1) independently of the xf ^ ~ !)• On^ the right hand 

side is an upsilon random variable. To perform the hypothesis test at the a 
level, compute the a quantile of the upsilon random variable with coefficient 



®And coincides with the lambda prime distribution. 
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and degrees of freedom v = [ni — l,n 2 — 1]^. If this quantile is bigger than 
zero, reject Hq in favor of Hi. 

Note that this two sample test can be further generalized to a fc independent 
sample test for a single equation involving signal-noise ratios. For i = 1,2,..., k, 
given Hi independent draws from Gaussian returns from k assets with signal- 
noise ratios Q, to test, for fixed ai, 02 ,..., a^, 6, to test the hypothesis 


Hn 


E 


diCi = b versus Hi 


E 


^ b. 


compute the sample Sharpe ratios, Q. As noted previously. 


c* 




and so, under the null. 




-1/2 





As previously, this random variable is an upsilon under the null, and so one re¬ 
jects at the a level if the left hand side exceeds the a quantile of the appropriate 
upsilon distribution. Note that this formulation subsumes the fc = 1 case. 

The analagous tests for the case of dependent returns does not seem to be 
simply defined. Rather one should rely on the standard asymptotic computation 
of the Sharpe ratio based on the central limit theorem and delta method, then 
test based on asymptotic normality, as described by Leung and Wong. [12, 19, 
18] 

The fc independent sample test can be applied to the case of Sharpe ratio 
under factor models as well.^ For i = 1, 2,..., fc, given ni independent draws 
from factor models with Gaussian errors on fc assets, for fixed Vi,V 2 , ■ ■ ■ ,Vk,c, 
to test the hypothesis 


Vi 

Hq : 2 ^ -= c versus 

O'i 


u Vi 

Hi : > - > c, 

^ ^ rr- 


compute the sample factor model Sharpe ratios for each sample: 


■5S,* 



Then check the probability that an upsilon distribution with coefficient 

Again it must be stressed that the factor returns should be deterministic, as otherwise 
their variability would contribute to extra uncertainty in the test statistics. 


-1 


- 1/2 




1 T 
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and degrees of freedom u = [ni — li,n 2 — h, - ■ ■ ,nk — h] exceeds 

- 1/2 




3.2 Confidence and Prediction intervals 


The one and two sample hypothesis tests can be rearranged to form confidence 
and prediction intervals. For example, a 1 — a confidence interval on C has 
endpoints [(^i , Cu] defined by quantiles of the upsilon distribution; 


VnCi = T„/2 

y/nCu = Ti_« 


( VnC , [n - 1]) , 

/2 ( VnC , [n - 1]) • 


( 11 ) 


By the same argument, an a/2 confidence level for the factor model signal-noise 
ratio, Cg, is given by [Cg,i,Cs,ti] defined as 


(«T(fTf)-^) '/^Cg/=T„/2(^(t;^(F^F)'^) Cg, [n - Z]) , 

1^2 ^ 1^2 '\ 
Ca,u = T,_^/y(v^{f^fy\y yin-y. 


Note that the confidence intervals use the k = 1 case of the upsilon, correspond¬ 
ing to the lambda prime distribution. 

Suppose, based on a sample of size ni, you observed Ci for some asset stream. 
What can you expect of the Sharpe ratio for n 2 future observations? The Fre- 
quentist answers this question via a prediction interval, an interval which, condi¬ 
tional on (/i and ni, contains the Sharpe ratio of those future observations with 
some specified probability, under replication.® Prediction intervals are actually 
an application of the two independent sample hypothesis test of equality. 

Suppose you observe Ci on ni observations of normally distributed i.i.d. 
returns, then observe (2 on n 2 observations from the same returns stream. As 
in the case of the two independent samples test, write 

Ci\lxl/{ni -1) + = C = C2\Jxl/ {^2 - 1) (13) 

where the Zi ^ Af {0, 1), and the Xi (^i ~ 1) are independent. Then, the 
probability that C 2 is less than some value, say ?/, is the probability that an 
upsilon with coefficient 


t = 


nin2 
ni -b 712 L 


Cl,-2/ 


and degrees of freedom zv = [ni — 1, n 2 — 1]^, is less than zero. 

^Typically ‘prediction interval’ is reserved for an interval around a single future observation, 
while ‘tolerance interval’ is used for multiple future observations. Our application is somewhat 
between these two. 


7 













Thus the prediction intervals are given by [C 2 ,i, C 2 ,tt]) defined implicitly by 


0 = T„/2 
0 = Ti_„/2 


nin2 
ni + 712 L 


1 T 


Cl,-C2.i ,[ni- 1,712-1] , 


711712 


ni + 712 1 


1 T 


Cij~C2,m ,[771 — 1,712 — 1] 


(14) 


To find these prediction intervals, one must resort to a numerical root finder 
applied to the cumulative distribution function of the upsilon. It should be 
noted that a more direct computation could be performed if one had access to 
the quantile function of Lecoutre’s K-prime distribution. [15, 21] Manipulating 
Equation 13, one has 


CivWK^ry+yWf 

Vxi/(7l2 - 1) 


(15) 


where Z is a standard normal independent of the two chi-squares. In 
Lecoutre’s terminology, conditional on observing ^i, (^2 is distribnted as 


K' 

ni — l,n 2 — 


1 (6,-/ 


+n2 ^ 


ni 
711712 


4 Bayesian inference on the Sharpe ratio 


The Frequentist tests considered previously all have similar Bayesian counter¬ 
parts. In the traditional development of Bayesian inference on a Gaussian dis¬ 
tribution with unknown parameters, prior and posterior distributions are con¬ 
sidered on the mean and variance, /r and cr^, or the mean and precision, the 
latter defined as cr“^. [10, 20] It is a simple task to reformulate these in terms 
of the signal-noise ratio, and some transform of, say, the variance. 

One commonly used conjugate prior is the ‘Normal-Inverse-Gamma’, under 
which one has an unconditional inverse gamma prior distribution on cr^ (this is, 
up to scaling, one over a chi-square), and, conditional on a, a normal prior on 
p.. [10, 3.3] These can be stated as 

cr^ oc T”^ (777 o/2, mcialli) , 

[i |(T^ CX M {po, cr'^/no) , 


where al,mQ,pQ and uq are the hyper-parameters. Under this formulation, an 
noninformative prior corresponds to mo = 0 = tiq. 

After observing n i.i.d. draws from a normal distribution, Af {p,a), say 
xi,X 2 , ■ ■ ■ ,Xn, let p and a be the sample estimates from Equation 4. The 
posterior is then 

cr^ oc T"^ (777 i/2, mial/2) , 

p\a‘^ (X Af [pi, (7^ /tii) , 


where 


77l = 77o -b 77, 


Ml 




nopo+np 
77l 

TTioo-o -b (77 - 1 ) 0-2 -b ^ {po - pf 

mi 


(18) 


mi = mo -b 77, 


(19) 



















This commonly used model can be trivially modified to one on the variance 
and the signal-noise ratio, where the former is a nuisance parameter. Trans¬ 
forming Equation 16, we arrive at 

cr^ oc T”^ (mo/2, mocro/2) , 

C|cr^ oc A/'(^^,l/no) , 

Marginalizing out cr^, we arrive at a upsilon prior 

y/n^C, oc T (v^Co, ?7io), 

where (/o = /^o/co- In this case the upsilon coincides with the lambda prime 
distribution. The marginal posterior can be written as 


( 20 ) 


( 21 ) 


V^C oc T (VnTCiinii): 


( 22 ) 


where 


n-i = no -b n, 


mi = mo -b n. 


Cl 


noCogo + nCd- 

niCTi 

moCTo + (n - 1) -b ^ (CoCo - 

mi 


(23) 

(24) 


where C, = jllu. 

The unattributed model can be generalized to the factor model by following 
the standard Bayesian regression analysis. Again, we are assuming that the 
factors are deterministic. The Bayesian regression prior is typically stated as 


cr^ oc r ^ (mo/2, mocrQ/2) , 
/3|cr^ oc A/'(/3 o,o'^Ao"^) , 


where , mo are the Bayesian hyperparameters for the coefficient and degrees 
of freedom of the error term, while /3 q is that for the regression coefficient, and 
Aq parametrizes uncertainty in the regression coefficient. 

Again assume one has n observations of stacked row-wise in the n x I 
matrix, F, and the corresponding returns stacked in vector x. Again define 

/3 =df (F''"F) ^F^£c and a =df \j— F/3^ {x — F/3^ (n — l)~^. The posterior 
distribution is 

cr^ oc T”^ {mil2, mialji) , 
ocAr(,3i,a2Ar^) , 

where 

Ai=Ao + FTf, /3i = Ai-i (Ao/3o + F^F/3) , 

2 moCT^-b (n - 0 +/3^F^F/3-b/3o^Ao/3o -/3 i^Ai/3i 

mi = mo + 71, ai = -. 

mi 

(28) 


(26) 

(27) 
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A non-informative prior corresponds to Aq = 0, = 0 , Ug = 0, mg = 0. 

We can collapse the prior or posterior ‘along’ the direction v via 

oc r-i (m,/2, , 

(3 \a^ (X Af Ai~^v) , 

where i = 0 for the prior and i = 1 for the posterior. As in the unattributed 
model, marginalizing out we have a upsilon prior and posterior: 


Ki ^v) = (t)^Aj ^v) (g (xT (^{v^Ai ^v) 


where Cg.i =df fBjvjai. 


(30) 


4.1 Credible and posterior prediction intervals 

One can construct credible intervals on the signal-noise ratio based on the pos¬ 
terior, so-called posterior intervals, via quantiles of the upsilon distribution. [10] 
For example, a (1 — a) credible interval on C is given by 

['l'a/2 (\AiiCi, mi), Ti_„/ 2 m-i)] , (31) 

where {t,u) is the q quantile of the upsilon distribution with coefficient t 
and u degrees of freedom. For the case of noninformative priors (corresponding 
to no = mg = 0), this becomes 




,T 


l-a/2 



(32) 


This is equivalent to the Frequentist confidence intervals given in 3.2, but re¬ 
placing n for n — 1 in the degrees of freedom for a. Asymptotically, the Bayesian 
credible interval for an noninformed prior is the same as the Frequentist confi¬ 
dence interval. Alternatively, a Bayesian quant could argue that her Frequentist 
cousin is a confused Bayesian with prior no = 0, mo = —1, = 0. 

For the factor model signal-noise ratio, the (1 — a) credible interval on C,g is 


1 

\/ Aiv 




i,mi 


, Ti _„/2 (^Vn^AivCg,! J 


(33) 


A Bayesian prediction interval is an interval which contains some fixed pro¬ 
portion of our posterior belief about the Sharpe ratio of some future observa¬ 
tions. This is very similar to the Frequentist prediction interval, but dances 
around the issues of frequency and belief that separate the Frequentist and 
Bayesian. 

As in the Frequentist case, our belief is that (' 2 , based on 712 future observa¬ 
tions, will be drawn from a compound non-central Sharpe ratio distribution with 
non-centrality parameter drawn from the posterior distribution. Effectively this 
is a ‘t of lambda prime’ distribution, as was the case in the Frequentist setting. 
We can summarize this as 


V^C oc T (^/nTCi, mi), 
\/^C 2 IC oc t (V^C, n 2 -l), 


(34) 
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although this jumbles up the usual notation, since C 2 is a quantity one can 
eventually observe, not a population parameter. Nevertheless, the intent of 
these equations should be clear. The consequence is that we can find posterior 
prediction intervals, as in the Frequentist case, by using the upsilon distribution. 
A 1 — a prediction interval on C 2 is given by [OcCft.*] where (lo is chosen such 
that the upsilon distribution with coefficient 


t = 


nin2 

Til + n2 


[Cl: -Oo]^, 


and degrees of freedom v = [mi,n 2 — 1]^, is less than zero with probability 
a/2, and C,hi is defined mutatis mutandis. 


5 An example: Momentum and the January ef¬ 
fect 

The methods above are illustrated using real asset returns, ignoring the as¬ 
sumption of Gaussian returns and deterministic factors. The monthly returns 
of ‘the Market’ portfolio, the small cap portfolio (known as SMB, for ‘small 
minus big’), the value portfolio (known as HML, for ‘high minus low’), and the 
momentum portfolio (known as UMD, for ‘up minus down’), as tabulated by 
Kenneth French, were downloaded from Quandl. [9, 24] These are the cele¬ 
brated three factor portfolios of Fama and French plus Car hart’s momentum 
factor. [8, 4] 

The data are distributed as monthly relative returns, quoted in percents. The 
Market return is quoted as an excess return, with the risk free rate subtracted 
out. The risk free rate is also tabulated. For our purposes, the raw market 
returns are needed, so the risk free rate is added back to the market returns. 
The set consists of 1056mo. of data, from Jan 1927 through Dec 2014. 

Now we will consider the hypothesis that the momentum factor has smaller 
Sharpe ratio in January than the remainder of the year, under a factor model 
where we attribute returns to the remaining three factors, a so-called ‘January 
effect’. We are effectively testing the information coefficient of UMD. Since the 
returns are independent®, we can apply the independent two sample test. 

Under the factor model, the Sharpe ratio for the Januaries was computed to 
be —0.137mo.“^/^, while for non-Januaries, it was computed to be 0.299mo.“^/^. 
To test the null hypothesis that the factor model signal-noise ratios are equal, we 
consider quantiles of the upsilon distribution with coefficient [—0.993, —2.175]^ 
and degrees of freedom [84,964]^. The 6 term Cornish Fisher approximate 
0.005 and 0.995 quantiles were computed to be [—5.751, —0.578], which does not 
contain zero. The 8 term Edgeworth approximation to the CDF was computed 
at zero to be 0.999, larger than 1 — a/2 with a = 0.01. Thus we reject the null 
hypothesis of equality at the 0.01 level in favor of the alternative hypothesis 
that the UMD portfolio has smaller information coefficient in January than the 
remainder of the year. 

®Independence ignores the small ‘bounce’ effect caused by any mispricing of assets in the 
portfolio at the time boundaries. This is likely to be small compared to the error in the other 
assumptions. 
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